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We investigate the periodic Anderson model with attractive interactions by means of dy- 
namical mean- field theory (DMFT). Using a continuous-time quantum Monte Carlo impurity 
solver, we study the competition between the superfluid state and the paramagnetic Kondo 
insulating state, and determine the phase diagram. At the chemical potential-induced phase 
transition from the Kondo insulating state to the superfluid state, a low-energy peak char- 
acteristic of the superfluid state appears inside the hybridization gap. We also address the 
effect of the confining potential in optical lattice systems by means of real-space DMFT 
calculations. 
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1. Introduction 

Ultracold atomic gases have attracted considerable interest 1-3 ) since the successful realiza- 
tion of Bose-Einstein condensation in a bosonic 87 Rb system. 4 ) One of the most active topics 
in this field is the study of fermionic optical lattice systems, which are formed by loading ultra- 
cold fermions into a periodic potential. 5-8 ) This setup provides a clean realization of a quantum 
lattice system, in which remarkable phenomena have been observed such as the Mott insu- 
lating state 9,10 ) and the BCS-BEC crossover in two component fermionic systems. 11-13 ) Due 
to the high controllability of the lattice structure and onsite interactions, ultracold fermionic 
systems in optical lattices can be regarded as quantum simulators of the Hubbard model, 
whose ground-state properties have been investigated by numerous theoretical and compu- 
tational approaches. Recent papers have suggested the possible realization of optical lattice 
systems described by other theoretical models such as multi-component models 14-16 ) and the 
Kondo lattice model. 17-19 ) This stimulates further investigations on lattice models with strong 
correlations. Among them, the periodic Anderson model (PAM) describing conduction and 
localized bands is one of the most important models in condensed matter physics. It captures 
the essential physics of some heavy-electron systems realized in rare-earth compounds. 20 ) An 
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important point is that quantum critical behavior is expected in the PAM with both repulsive 
and attractive interactions, and the model may thus provide a stage to discuss a quantum 
phase transition in an optical lattice system. However, this type of instability, which may 
be important to predict the low-temperature properties of optical lattice systems, has not 
been discussed in the PAM with attractive interactions. Furthermore, the confining potential 
should play a crucial role in understanding the quantum critical behavior in optical lattice 
systems with attractive interactions. Therefore, it is important to systematically investigate 
the low temperature properties in this model with localized and itinerant bands. 

In this paper, we study the PAM with attractive interactions by combining dynamical 
mean-field theory (DMFT) 21 ~ 24 ) with the continuous-time quantum Monte Carlo (CTQMC) 
approach. 25 ) First, we discuss how the superfluid state is realized on a uniform bipartite lattice 
and determine the phase diagram. By examining the dynamical properties in the superfluid 
state in detail, we demonstrate that low energy in-gap states are induced in the density of 
states. We also discuss the effect of the confining potential by means of the real-space DMFT 
and investigate how the Kondo insulating state spatially competes with the superfluid state. 

The paper is organized as follows. In Sec. 2, we introduce the model Hamiltonian and 
explain the particle- hole transformation in the PAM. We briefly summarize our theoretical 
approach. We demonstrate how the superfluid state competes with the Kondo insulating state 
at low temperatures in Sec. 3. The effect of the confining potential is discussed in Sec. 4. A 
brief summary is given in the last section. 

2. Model and Method 



We consider the periodic Anderson model with conduction and localized bands, which is 
described by the following Hamiltonian, 



H 
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where Ci a (fi a ) annihilates a particle in the conduction (localized) band on the ith site with 
spin cr, and nf a = a.\ (a = c and /). t is the nearest-neighbor hopping matrix for the 
conduction band, V is the hybridization between the two bands, U is the attractive interaction, 
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and Ef is the energy level for the localized band. fj, c (fj,f) and h c (hf) are the chemical potential 
and the magnetic field for the conduction (localized) band. In the following, we refer to the 
conduction and localized bands as c and / bands, for simplicity. 

First, we wish to note that the low-energy properties of the PAM with attractive interac- 
tions [eq. (1)] are closely related to those of the repulsive PAM. This can be understood by 
applying a particle-hole transformation to the model on the bipartite lattice. It is explicitly 
given by 



c ia -> ac\ 



i a 



i g A, (4) 



3 e B, (5) 



fj* -»• a fja 

where A and B are sublattice indices, and Ci a (fi a ) annihilates a particle in the c (/) band on 
the ith site with spin a. We then obtain the PAM with repulsive interactions U(= —U), where 
the energy level for the / band, chemical potentials and magnetic fields are transformed as 
Ef = 0, [if = hf, jjj C = h c , hf = Ef — [if, and h c = fi c , respectively. In particular, the low 
temperature properties of both models are equivalent at half filling (Ef = 0, fj, a = 0, h a = 0), 
which is similar to the case of the Hubbard model. 26 ) 

The low-temperature properties of the PAM with repulsive interactions, which have been 
discussed in the context of / electron systems, 20 ) may be helpful to understand those of our 
attractive model. In the repulsive PAM, the hybridization between the c and / bands favors 
the formation of a local Kondo singlet state. Therefore, a large hybridization stabilizes the 
Kondo insulating state at half filling and the heavy metallic state away from half filing. 27 ' 28 ) 
On the other hand, the Coulomb repulsion favors magnetic correlations between sites, which 
tends to induce an antiferromagnetically ordered state. 29 ~ 33 ) Therefore, in the half-filled PAM 
on a bipartite lattice, a quantum phase transition occurs between the Kondo insulating state 
and the antiferromagnetically ordered state in two and higher dimensions. Note that this 
competition is controlled by not only the interaction strength but also the magnetic field. If 
a magnetic field is applied, the spin gap in the Kondo insulating state decreases and finally 
vanishes. At this point, a field-induced magnetic phase transition occurs, where a spontaneous 
staggered magnetization appears in the plane perpendicular to the applied field. 34-36 ) 

In the PAM with attractive interactions, a similar competition is expected at low tem- 
peratures. The corresponding ordered phase should be characterized by a pair potential A, 
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which is given by 

Aa = -^^2( a n a n)' (6) 

i 

where a = c, /, and N is the total number of sites. This quantity corresponds to the staggered 
magnetization in the x direction = + ^L^if) m the PAM with repulsive 

interactions. By contrast, it is known that the Kondo insulating state has no order parameter 
since it is adiabatically connected to the paramagnetic state at high temperatures. This implies 
that it is also realized in the PAM with attractive interactions. 

To discuss how the superfluid state competes with the Kondo insulating state, we make use 
of DMFT. Since local particle correlations can be taken into account precisely, this method is 
formally exact in the limit of infinite dimensions, and has successfully been applied to strongly 
correlated systems. In DMFT, the lattice Green's function is obtained via a self-consistency 
condition imposed on the impurity problem. The non-interacting Green's function for the 
lattice model is given as, 

(V (MO = (7) 
\ -Vcr z £/er + (pf - E f ) a z J 

where £ a = iuj n +h a (a = c, /), is the dispersion relation of the c band, a z is the z component 

of the Pauli matrix, <ro is the identity matrix, oj n = (2n + 1)ttT is the Matsubara frequency, 

and T is the temperature. Here, the Green's function is represented in the Nambu formalism 

to describe the superfluid state. The lattice Green's function is then given in terms of the site 

diagonal self-energy T,(iuj n ) as, 

G c (iu! n ) G c f(iu) n ) 
G(ico n ) = | cK n> /V ^ | (8) 
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S(iw n ) = ( | . (9) 

Ef(iu) n ) 

The self-energy for the / band Sj is given as, 

m f\ i (10) 

Sf{iu n ) -E* f ^(iu n ) 

where Y,f a (iuj n ) [Sf(iui n )] is the normal (anomalous) part of the self-energy for the / band. 
In DMFT, the self-consistency condition is given by 

G f (iu n ) = G imp (iuj n ), (11) 
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where Gi mp is the Green's function of the effective impurity model. The effective medium for 
each site is given by 

g-H^n) = [Gf (iiOn)}' 1 + £ (iu n ) . (12) 

In the PAM with attractive interactions, a competition between the Kondo insulating 
state and the superfluid state is expected. Therefore, it is necessary that the impurity solver 
is capable of treating different energy scales accurately. The exact diagonalization method 37 ) 
is efficient to discuss ground state properties, but it is difficult to quantitatively discuss the 
critical behavior and dynamical properties at finite temperatures. The numerical renormaliza- 
tion group method 38, 39 ) is one of the most powerful methods to describe the low energy states, 
but may not be appropriate for investigating the higher energy region. In this paper, we use 
the recently developed CTQMC method. In this approach, a Monte Carlo sampling of certain 
collections of diagrams is performed in continuous time, and thereby the Trotter error, which 
originates from the Suzuki- Trotter decomposition, is avoided. The CTQMC method comes in 
two flavors, a weak-coupling 25 ) and strong-coupling 40 ) formalism, and is applicable to general 
classes of models such as the Hubbard model, 41 ~ 46 ) the periodic Anderson model, 47 ) the Kondo 
lattice model, 48 ) and the Holstein-Hubbard model. 49 ) Here, we use the continuous-time auxil- 
iary field version of the weak-coupling CTQMC method 50 ) extended to the Nambu formalism, 
which allows us to directly access the superfluid state at low temperatures. 45 ) In our CTQMC 
simulations, we measure normal and anomalous Green's functions on a grid of a thousand 
points. To discuss static and dynamical properties in the system quantitatively, the number 
of CTQMC samplings in each DMFT iteration is changed, depending on the temperature, 
the chemical potential, and the magnitude of the hybridization: e.g. we have performed eighty 
million (five billion) samplings for the half-filled system with U = 1, V = 0.40, and T = 0.02 
(U = 1, V = 0.18 and T = 0.01). In the following, we set /i = [i c = /i/ and h = h c = hf = 0, 
for simplicity. 

3. Periodic Anderson model on the hypercubic lattice 

We first consider the PAM on the infinite-dimensional hypercubic lattice. The bare density 
of states for the c band is Gaussian with the bandwidth t*: p(uj) = exp[— (uj/t*) 2 ]/V irt* 2 . In 
this section, we use t* as the energy unit. The normal state properties have been discussed 
in detail by means of DMFT, and the phase diagram has been obtained from the divergence 
of the staggered susceptibility at finite temperatures. 30 ~ 33 ) However, the off-diagonal ordered 
state has not been discussed directly since it cannot easily be treated. Here, we combine DMFT 
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with the CTQMC method based on the Nambu formalism to discuss how the superfluid state 
is realized at low temperatures. In Fig. 1, we show the temperature dependence of the pair 
potential in the half-filled system with U = 1 and V = 0.3. At high temperatures, the pair 




T 

Fig. 1. (Color online) The pair potentials A c and Af for the superfluid state as a function of tem- 
perature T when U = 1 and V = 0.3. The inset shows the critical behavior of the pair potentials. 
Solid lines are guides to eyes. 



potential is zero and the normal state is realized. As temperature is decreased below a critical 
value, the pair potentials for both bands are simultaneously induced, which implies a phase 
transition to the superfluid state. By examining the critical behavior A ~ \T — T c \@ with the 
exponent f3 = 1/2, we determine the critical temperature T c ~ 0.023, as shown in the inset of 
Fig. 1. 

However, the superfluid state is not always realized at low temperatures in the PAM. By 
performing similar calculations for several values of the hybridization V, we find, as shown in 
Fig. 2, three distinct regions: V < V c \, V c \ < V < V C 2 and V C 2 < V. The two critical points 
are obtained as V c \ ~ 0.26(0.18) and V c2 ~ 0.36(0.37) at T = 0.02(0.01). To clarify the nature 
of these states, we show the density of states obtained by the maximum entropy method in 
Fig. 3. When V > V C 2, the pair potential is zero and a charge gap appears around the Fermi 
level in each band. This implies that the large hybridization stabilizes the Kondo insulating 
state at half filling. In the intermediate region (V c i < V < V C 2), the superfluid state is realized 
since the system has finite pair potentials and the corresponding gap in the density of states. 
Note that no drastic change in the density of states appears at V = V C 2, in contrast to the 
pair potential. A similar behavior has been discussed in the half-filled Kondo lattice model. 51 ) 
On the other hand, when V < V&, a finite density of states appears around the Fermi level 
in each band. This implies that a normal metallic state is realized. 
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Fig. 2. (Color online) The pair potential as a function of the hybridization V when T = 0.01 (open 
symbols) and T = 0.02 (solid symbols). 



By performing similar calculations, we obtain the phase diagram for the half filled system 
with U = 1, as shown in Fig. 4. The normal state is realized at higher temperatures, while 
a different behavior appears at lower temperatures, depending on the magnitude of the hy- 
bridization. In the small V case, the system is reduced to the non-interacting c band weakly 
coupled to the correlated / band, which makes the critical temperature for the superfluid 
state rather low. Therefore, the normal metallic region shrinks as temperature is decreased, 
V& — > (T 0), as shown in Fig. 4. On the other hand, we find that the other critical 
point V C 2 between the superfluid state and the Kondo insulating state is slightly increased as 
temperature is lowered. We expect that at zero temperature, the superfluid state is realized 
in the small V region and a quantum phase transition occurs to the Kondo insulating state 
at V^2(~ 0.38). This is consistent with ground-state properties for the PAM with repulsive 
interactions 31 ) and the Kondo lattice model, 42 ' 48 ) where the Kondo insulator competes with 
the magnetically ordered state. 

Now, let us consider how the Kondo insulating state competes with the superfluid state 
away from half filling. In the following, we show only the pair potential A f since A c behaves es- 
sentially the same. In Fig. 5, we show the local particle density at each site n(= 'Ylaia^io) /-^O 
and the pair potential Af when the chemical potential is varied. Here, we show only the local 
particle density at T = 0.02 in Figs. 5 (a), (b) and (c) since n(T = 0.01) is almost identical. In 
addition, the compressibility k(= dn/dfj) deduced from the numerical derivative is also shown 
in Fig. 5 (d). This quantity corresponds to the magnetic susceptibility in the z direction in 
the repulsive PAM. Therefore, a cusp singularity is naively expected at the critical point. 
When V = 0.3, the half-filled system (/x = 0) is in the superfluid state at low temperatures, 
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Fig. 3. (Color online) Density of states for localized and conduction bands when V = 
0.25, 0.3, 0.35, 0.4 and 0.5 at T = 0.02. 



as discussed above. Shifting the chemical potential away from zero, the pair potential is little 
affected and the local particle density varies smoothly, as shown in Fig. 5 (a). Therefore, we 
conclude that the superfluid state is stable against small changes in the chemical potential. 
Far away from half filling, the pair potential decreases and finally reaches zero. A phase tran- 
sition then occurs to the normal metallic state. By contrast, a different behavior appears in 
the case V = 0.35. The pair potential initially increases as the chemical potential is shifted 
away from zero, has a maximum value around (i ~ 0.05, and finally vanishes at the phase 
transition, as shown in Fig. 5 (b). The nonmonotonic behavior is even more clearly evident 
in the case V = 0.4, as shown in Fig. 5 (c). When fj, < fj, c \, where [i c \ ~ 0.049(0.032) at 
T = 0.02(0.01), the pair potential is zero and the local particle density is little changed. In 
fact, Fig. 5 (d) shows the strong suppression in the compressibility. This implies that an in- 
compressible Kondo insulating state is realized at half filling, instead of the superfluid state. 
An increase in the chemical potential beyond \x c \ drives a chemical potential-induced super- 
fluid phase transition, marked by the appearance of the pair potential and a cusp singularity 
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Fig. 4. (Color online) Phase diagram of the half-filled PAM with U = 1. Crosses are the critical points, 
solid (open) circles indicate the state with (without) the pair potential. The phase boundary is a 
guide to the eyes. 




Fig. 5. (Color online) The local particle density and the pair potential for the / band as a function of 
the chemical potential when V = 0.3 (a), V = 0.35 (b), and V = 0.4 (c). (d) The compressibility 
of the system with V = 0.3, 0.35 and 0.4 at T = 0.02. 



in the compressibility although the latter is not clearly visible in Fig. 5 (d). Further increase 
in \x changes the local particle density and finally drives the system to the normal metallic 
state at // = fi c2 , where (i c2 ~ 0.10(0.14) at T = 0.02(0.01). 

By performing similar calculations for various temperatures, we obtain the phase diagram 
for U = 1 and V = 0.4 shown in Fig. 6. This is essentially the same as the phase diagram for 
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Fig. 6. (Color online) Phase diagram for the system with U = 1 and V = 0.4. Crosses are the 
critical points, solid (open) circles indicate the state with (without) the pair potential. The phase 
boundary is a guide to the eyes. 



the magnetic field in the PAM with repulsive interactions. 34 ~ 36 ) When the system is half filled 
(fi = 0), the Kondo insulating state is realized at low temperatures. At T = 0.01, two peak 
structures clearly appear at the edges of the hybridization gap in the density of states of the 
/ band, as shown in Fig. 7 (a). The introduction of the chemical potential simply shifts these 



(a) Pf(co) 




Fig. 7. Density of states for the / (a) and c (b) bands of the system with U = 1 and V = 0.4 at 
T = 0.01. The data are for the chemical potentials fx = 0.00,0.01, • • ■ ,0.17 from the bottom to 
the top. 



peaks in the case \i < 0.03 since the Kondo state is incompressible. Further increase leads to 
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interesting behavior in the density of states. We find that around fi ~ 0.05, the upper peak 
never approaches the Fermi level and low energy states are, instead, induced below the Fermi 
level between the two peaks. This implies that particles in the vicinity of the Fermi level form 
Cooper pairs, which yields low-energy in-gap states. The latter interpretation is consistent with 
the fact that the chemical potential-induced superfluid phase transition occurs at jjL c \ ~ 0.032. 
When the chemical potential increases beyond the critical point, the low-energy feature grows 
and becomes the dominant peak. It is also found that as the chemical potential is varied, the 
center of the gap is always located at the Fermi level, in contrast to the case with [i < [i c \. 
On the other hand, the peak at the lower edge of the hybridization gap becomes smeared. 
Therefore, we conclude that a compressible superfluid state is stabilized in this region. When 
the system approaches /x C 2(~ 0.14), the superfluid gap collapses and another phase transition 
occurs to the normal metallic state. A similar behavior is also found in the density of states 
of the c band, as shown in Fig. 7 (b). We emphasize that our calculation directly treats the 
low temperature superfluid state, in contrast to previous works for the PAM. 34 ~ 36 ) Figure 7 
therefore clarifies how the chemical potential-induced phase transition yields low energy in-gap 
states in the density of states. 

The results obtained above may have important consequences for three-dimensional optical 
lattice systems. In particular, one might expect quantum critical behavior in the vicinity of 
a certain surface in the three-dimensional optical lattice system if the confining potential is 
regarded as a site-dependent chemical potential. In the following section, we study the effect 
of the confining potential carefully to discuss the competition between the Kondo insulating 
state and the superfluid state in such optical lattice systems. 

4. Effect of the confining potential 

We study here the effect of the confining potential in the PAM. As discussed in the previous 
section, we have found that the Kondo insulating state competes with the superfluid state 
at low temperatures. Furthermore, it has been clarified that in a certain parameter region, a 
shift of the chemical potential can stabilize the superfluid state. Therefore, if the confining 
potential is simply regarded as a site-dependent chemical potential, the Kondo insulating 
state and the superfluid state may be separated by a quantum critical surface in the three- 
dimensional optical lattice system. In this case, the critical temperature should depend on the 
position in the trap. However, the physics in the presence of a trap is highly nontrivial and 
therefore it is necessary to discuss the low temperature properties using more sophisticated 
tools. 
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To this end, we use the real-space DMFT. The method is based on the local approximation, 
but intersite correlations are to some extent taken into account. In fact, the real-space DMFT 
approach has been successfully used to discuss low temperature properties of optical lattice 
systems 46 ' 52 ~ 54 ) and the interface between transition metal oxides. 55 ) Here, making use of 
this method, we consider a three-dimensional optical lattice system. The Hamiltonian for the 
confining potential is defined as, 

H c = J>(ri)rC (13) 

iaa 

where v(r) = vo(r/a) 2 , vq is the curvature of the confining potential, r% is the distance between 
the ith site and the center of the trap, and a is the lattice constant. In this section, we fix 
t = 1 and vq = 0.1. We then treat a system with 4169 sites restricted by the condition 
{r/a < 10). Using the CTQMC method as an impurity solver, we perform the real-space 
DMFT calculations. We obtain results for the spatial distributions of the local particle density 
and the pair potentials for the system with U = 3.5 and fi = 0.5, as shown in Fig. 8. When the 




Fig. 8. (Color online) Local particle density (a) and the pair potentials (b) as a function of the 
distance from the center in the three-dimensional optical lattice system. Triangles, circles, and 
squares represent the results in the system with V = 1.0,1.5 and V = 2.0 when U — 3.5 and 
T = 0.033. The inset of figure (a) shows the local particle density for the c and / bands. The pair 
potential is shown on a logarithmic scale in the inset of the figure (b). 
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hybridization is small (V = 1.0), the localized / particles are weakly coupled to the c band. 
Therefore, the / band is almost occupied in the core r/a < yjtjvo ~ 2.2, and empty outside 
the core due to the attractive interaction for the / band. By contrast, the noninteracting c 
particles are widely distributed due to the hopping integral between sites, as shown in the 
inset of Fig. 8 (a). In this case, no pair potential appears and the metallic state is realized. 
An increase in the hybridization affects the distribution of the / band and low-temperature 
properties. When V = 1.5, the local particle density around the center is decreased and the 
pair potential is induced in the system, as shown in Fig. 8. Therefore, the superfiuid state is 
realized in the whole region with n(r) ^ although the pair potential is tiny away from the 
center, as shown in the inset of Fig. 8 (b). This behavior is similar to that in the attractive 
Hubbard model with a confining potential. 46 ^ Note that the rapid decay in the pair potential 
sets in around r/a = 1, in contrast to the distribution of the local particle density. This may 
be explained by the competition between the Kondo insulating state and the superfiuid state. 
We find a shoulder structure around r/a = 2 in the local particle density, which indicates a 
tendency to form the Kondo insulating state in the superfiuid state. This results in the rapid 
decay of the pair potential. Further increase in the hybridization tends to result in a plateau 
structure and decreases the pair potential. When V = 2.0, the pair potential vanishes and 
the plateau clearly appears at a value of n(r) corresponding to half of the states occupied. 
Therefore, the Kondo insulating state is stabilized in the core region with r/a< 3. 

With these calculations, we have clarified that the Kondo insulating state spatially com- 
petes with the superfiuid state in an optical lattice system with confining potential. However, 
we could not find a critical surface between the superfiuid state and the Kondo insulating state, 
which could naively be expected from the results of the infinite dimensional PAM. Therefore, 
we can say that intersite correlations play an important role and must be considered in a 
discussion of the low temperature properties of optical lattice systems. 

Before concluding the paper, we would like to comment on correlation effects in the c 
band, which can be tuned by the Feshbach resonance in optical lattice systems. According to 
previous studies of the PAM with repulsive interactions, these interactions do not affect the 
low temperature properties qualitatively. 56,57 ) Therefore, it is expected that in the attractive 
case, the essence of the low temperature properties must be described by the simple PAM [eq. 
(1)] although phase boundaries may be slightly shifted by interactions between the c particles. 
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5. Summary 

We have investigated the periodic Anderson model with attractive interactions on the 
hypercubic lattice. By combining DMFT with the CTQMC method based on the Nambu 
formalism, we have studied quantitatively how the superfluid state is stabilized at low tem- 
peratures. It has been found that a low-energy state characteristic of the superfluid phase 
appears in the hybridization gap near the critical chemical potential. We have also discussed 
the effect of the confining potential in the three-dimensional optical lattice by means of the 
real-space DMFT to clarify how the superfluid state spatially competes with the Kondo insu- 
lating state. 
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